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The spatial configuration of initial partons in high-multiplicity proton-proton scatterings at 14 
TeV is assumed as three randomly positioned "hot spots". The parton momentum distribution in 
the hot spots is calculated by HIJING2.0 with some modifications. This initial condition causes not 
only large eccentricity £2 but also triangularity £3 and the correlation of £2 — £3 event-plane angles. 
The final elliptic flow V2, triangular flow 1)3, and the correlation of v% — V3 event-plane angles are 
calculated by using the parton cascade model BAMPS to simulate the space-time parton evolution. 
Our results show that the V2 — 1)3 correlation is different from that of £2 — £3 • This finding indicates 
that translations of different Fourier components of the initial spatial asymmetry to the final flow 
components are not independent. A dynamical correlation between the elliptic and triangular flow 
appears during the collective expansion. 

PACS numbers: 



I. INTRODUCTION 

A strong interest in high multiplicity events in ultra- 
relativistic proton-proton collisions has arisen recently, 
since a near side "ridge" has been found in such events 
by the CMS collaboration at the Large Hadron Collider 
(LHC) at y/s = 7 TeV [l|. This phenomenon is very sim- 
ilar to that observed in Au-Au collisions at the Relativis- 
tic Heavy Ion Collider (RHIC) at y/sj^ = 0.2 TeV 
Because the quark-gluon plasma (QGP) created at RHIC 
behaves like a nearly perfect fluid one may suppose 
that a similar hydrodynamic behaviour has appeared in 
high multiplicity p-p collisions at LHC @. Methods de- 
veloped to investigate QGP at RHIC can be used to pre- 
dict new phenomena in p-p collisions at LHC. In this 
work we are concentrating on the elliptic and triangular 
flow induced by the initial eccentricity and triangularity 
in ultrarelativistic high multiplicity p-p collisions at the 
highest LHC energy. 

The elliptic flow parameter V2 is the best experimental 
observable determining the strength of the hydrodynamic 
collectivity 0-Qll • I n the presence of strong interactions 
of system constituents, i>2 is obtained by the translation 
of the spatial asymmetry of initially produced matter into 
the final particle angular distribution [ll|, EH- There 
are two ways to make a spatial asymmetry in nucleus- 
nucleus collisions. One is the geometric overlap in non- 
central collisions. This has been adopted to p-p collisions 
at LHC, because a proton also has its extension though 
small. However, the predictions following the geometric 
overlap-eccentricity £2 showed small values of the elliptic 
flow parameter V2 ~ 3% in minimum bias, either in hy- 
drodynamic calculations [ili - fl5| or in the £2 — v 2 scaling 
The latter is assumed to be the same as that pro- 
posed for nucleus-nucleus collisions at RHIC flR [l9j . In 
the geometric picture non-central collisions provide large 



initial eccentricity, which is the necessary condition for 
large elliptic flow. However, the particle multiplicity in 
such collisions in the case of p-p scatterings is rather low. 
Nonflow effects can be hardly eliminated. 

The second source of a spatial asymmetry comes from 
statistical density fluctuations of the initially produced 
matter on the event-by-event basis [2(|. This is the rea- 
son for the nonvanishing V2 (and also higher harmonics) 
in central nucleus- nucleus collisions at LHC [2l[ . Due to 
the much smaller volume of a proton compared with a 
Au nucleus, it is natural to consider fluctuations in cen- 
tral p-p collisions at the LHC energy, which can provide 
both initial eccentricity and high multiplicity. A few hot 
spots or flux tubes may be excited in a p-p collision and 
lead to measurable elliptic flow at LHC |22h25| . More- 
over, parton evolutions and multiple scatterings [2H [27j 
can also generate large event- by- event fluctuations. 

Unlike the smooth initial distribution that generates 
only even-order Fourier components in the momentum 
angular distribution, initial fluctuations lead to nonva- 
nishing odd-order components, which are shown to con- 
tribute the azimuthal correlations observed at RHIC [28r - 
3lj . In particular, the triangular flow V3 from the initial 
triangularity £3 is intensively studied in the recent works 

32 - |38{ . In central Au+ Au collisions U3 is as large as i>2 

371. 



For p-p collisions at LHC we suggest that V3 is as im- 
portant as V2 in high multiplicity events. Different from 
Au+Au collisions at RHIC, where calculations 0, |33[ 
indicate no correlations between the initial event-plan an- 
gles of £2 and £3 and between the final event-plan angles 
of V2 and W3, such correlations could exist in the p-p col- 
lisions at LHC due to the smaller number of hot spots. 
Measurements on the flow correlations will shed light on 
details in the context of collective flow phenomena. In 
this work we calculate £2, £3, «2, ^3, and their correla- 



2 



tions. The event-by-event generation of the parton initial 
conditions is an implementation of the hot spots scenario 
[li] by using the recent version of HIJING2.0 (2§E[. 
The dynamical space-time evolution is calculated by us- 
ing the parton cascade BAMPS [42j|. 



II. INITIAL CONDITIONS 



HIJING [39H41j is a Monte-Carlo event generator 
for hadron productions in high energy nucleon-nucleon, 
nucleon- nucleus, and nucleus-nucleus collisions. It is es- 
sentially a two-component model, which describes the 
production of hard parton jets and the soft interaction 
between nucleon remnants. While the hard jets produc- 
tion can be calculated by perturbative QCD (pQCD), 
nucleon remnants interact via soft gluon exchanges de- 
scribed by the string model [43| . The produced hard jet 
pairs and the two excited remnants are treated as inde- 
pendent strings, which fragment to resonances that decay 
to final hadrons. The predictions using the updated HI- 
JING2.0 [Hi are in good agreements with the recently 
measured hadron spectra at LHC in p+p collisions at 
Vs = 0.9, 2.36 and 7 TcV 0-51, and central Pb+Pb 
collisions at ^/sjvjv = 2.76 TeV |49j . 

Using HIJING2.0 we calculate the hadron multiplic- 
ity for p+p collisions at 14 TeV and find that the total 
hadron multiplicity dN/dy at y = has a mean value 
of 10.4 and possesses a high- multiplicity tail reaching 
the abundance of semi-peripheral Cu+Cu collisions at 
t/sn~n = 62 .4 GeV at RHIC [5Cj. In this work we are in- 
terested in the events in the window of 50 < dN/dy(y — 
0) < 60. In these events the probability to produce three 
strings (one is formed by the hard jet pairs and other two 
from the excited remnants) is more than 80%. We thus 
neglect, for simplicity, events without hard jets (17%) 
and events containing more than one jet string (3%). 
The rapidity distributions of the hadron multiplicity and 
transverse energy in these selected events are shown in 
Fig. [TJ The peak at midrapidity is due to the fact that 
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FIG. 1: (Color online) The rapidity distribution of multiplic- 
ity and transverse energy of hadrons produced in high mul- 
tiplicity p+p collisions at 14 TeV, calculated by HIJING2.0 
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in the selected events, the hard jets production and the 
multiple gluons exchange in the soft interaction mostly 
occur at midrapidity. 

To generate the initial condition on the parton level, 
we turn off the resonance decays in HIJING and return to 
the representation of resonances by quark-antiquark pairs 
or q uark- diq uark pairs according to the LUND string 
breaking |43l . [5l| . This approach is similar to the string 
melting implemented in AMPT [52| , where final hadrons 
are converted into partons. Figure shows the rapid- 
ity distribution of parton number and transverse energy 
from jet, projectile, and target string, respectively. 
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FIG. 2: (Color online) The same as Fig. [TJ but on the parton 
level from three strings. See details in text. 

Because the three strings break independently, we as- 
sume that partons from each string build a hot spot. Due 
to the high collider energy all partons are produced at 
2 = 0. The spatial distribution of the three hot spots in 
the transverse plane follows the scenario proposed in Rcf. 
[23| . The center of the hot spots is determined according 
to the proton spatial density [53[ 
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where n = 0.17/ /m 3 , R = 0.56 fm and d = 0.112 
fm. The spatial parton distribution in each hot spot is 
assumed to have a Gaussian profile e~ r / r ° , where ro 
gives the size of the hot spots. 

Although r of the remnant could be larger than that 
of the jet string, we assume that r is equal for each hot 
spot, which gives the largest effect on the initial asym- 
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to be 100. Using exact numbers from HIJING (about 160 
in each hot spot taken over all rapidities, see Fig. 2) has 
tiny changes on the results. 100 partons per hot spot are 
roughly the numbers within the rapidity range \y\ < 2.5, 
which is a similar region covered in CMS [l| . Our results 
show that the most 63 arc not much smaller than the 
most 62, which indicates the possibility to measure both 
vi and U3 experimentally, if the assumptions of the initial 
condition and the hydrodynamic transport arc justified. 

Due to the random nature of the initial condition, the 
distribution of $2 and $3 are uniform within the interval 
of [— 7r/2,7r/2] and [— 7r/3, 7r/3], respectively. The lower 
panel of Fig. [3] shows the probability distribution of the 
angle difference, A$ = | <I> 2 — ^3 1 , which indicates the 
correlation of the two initial event-planes. Events with 
A$ = 30 degree are more probable. 

For initial conditions with two hot spots instead of 
three, we expect vanishing £3. With four or more hot 
spots both 63 and ti become smaller and their event- 
plane angles will be less correlated. Therefore, the largest 
correlation comes from events with three hot spots. 

When hot spots expand and overlap, collective flow 
components in higher order will be built up. By analogy 
with e„ flow coefficients v n are defined as the harmonic 
components of the particle azimuthal distribution in mo- 
mentum space 



FIG. 3: (Color online) Upper: The probability distribution 
of spatial eccentricity and triangularity. Lower: Correlation 
(arbitrary unit) of the e 2 — £3 event-plane angles. 



metry. is set to be 0.2 fm. The smaller the ro, the 
larger the e2 and 63 [23l |. 

Translating particles into the frame where the average 
position is equal zero, i.e., (x) = (y) = 0, the harmonic 
components e n of the spatial azimuthal distribution are 
defined as 133 
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where r, <fi are parton polar coordinates. The correspond- 
ing initial event-plane angles are given by 
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These are the angles $ = $„ where (— r™cosn(</> — 
$))/(r") has the maximum, which is e„ [see Eq. fl2J)]. 
For instance, for an ellipse shape with the short axis in 
the x direction <J> 2 = and £2 = (y 2 — x 2 )/(y 2 + x 2 ). 
Our definition of $„ differs from that from Ref. [33| by 
±ir/n. 

The upper panel of Fig. [3] shows the probabilities of 
the eccentricity e 2 an d the triangularity e 3 in the high- 
multiplicity p+p collisions at 14 TcV on the event-by- 
event basis. The parton number in each hot spot is set 
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is the event plane angle [54(, i.e., (cosn(ip 

= has the maximum, which is v n . 

If the translations from e„ to v„ for all components 
are completely independent, we will obtain "$! n = 
The strong correlation of A$ seen in Fig. [3] will be ob- 
served in A^ = I ^>2 — ^3] too. If such correlation can 
be measured experimentally at LHC, this will be the ev- 
idence for the hot spot scenario of the initial condition 
and the hydrodynamic behaviour of the parton matter in 
high-multiplicity p+p collisions. In this work we discuss 
this issue by simulating the parton collectivity within a 
microscopic manner. 



III. ELLIPTIC, TRIANGULAR FLOW AND 
THEIR CORRELATION 



The space-time evolution of the partons is simulated 
by the parton cascade model BAMPS (Boltzmann Ap- 
proach of MultiParton Scatterings) [42[ , which solves the 
Boltzmann equation for on-shell partons. For simplicity 
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we regard partons, which stem from the string break- 
ing (see Sec. [TTJ) , as identical massless Boltzmann par- 
ticles. The particle degeneracy factor is assumed to be 
the same as that of gluons together with quarks with 
two flavours. This leads to the particle number density 
n eq = (40/7r 2 )T 3 , where T is the temperature, if the sys- 
tem is in local thermal and chemical equilibrium. Fur- 
thermore, we consider only elastic binary scatterings and 
assume the isotropic distribution of the collision angle. 

The Boltzmann equation applies for systems when the 
particle mean free path A m / P = l/(na) is larger than the 
mean particle distance d = n -1 / 3 , where n is the local 
particle number density. For chosen constant \ m fp/d ra- 
tio as a global parameter we determine the total cross 
sections a in local cells, which are used to simulate scat- 
terings (42j . 

To assess hydrodynamic behaviour of the partons con- 
sidered in this work we calculate the shear viscosity to 
the entropy density ratio rj/s for given X m fp/d ratio. 
The shear viscosity is proportional to the energy density 
and the mean free path [55[, r\ = (2/5)eA TO / p . Assum- 
ing local thermal equilibrium, the entropy density is s = 
(4 — In A p )n, where A p = n/n eq denotes the particle fugac- 

1 /3 

ity. We have then r\js = 0.752A m / p /(i/(4 — lnA p )/A p . 
The relation e = 3nT is used in the last equation. The ra- 
tio rj/s has a weak dependence on the fugacity. Therefore, 
we take rj/s = 0.188A m / p /<i, which is exact for A p = 1, as 
an estimate of the rj/s value for systems out of equilib- 
rium, such as the present case. For choosing X m f p /d = 2 
for instance we obtain rj/s « 0.376. 

Further model parameters are set as follows: The ini- 
tial time for starting BAMPS is 0.1 fm/c. Before that 
time partons propagate freely. The cell length in the 
transverse plane is Ax = Ay = 0.02 fm, while the lon- 
gitudinal cell size is Arj w 0.1, expressed in the space- 
time rapidit y 77 = 0.51n[(t + z)/(t — z)]. The test parti- 
cle method [42j is used to enhance the numerical accu- 
racy. For that the parton density is amplified by a factor 
of 3000. Cross sections are reduced by the same factor 
to keep the mean free path unchanged. Parton scatter- 
ings stop when the local energy density is lower than 
lGeV/fm 3 . This mimics the phase transition, which has 
to be implemented in BAMPS in the future. For the 
event-by-event analysis we compute ten thousands runs. 

Figure [1] shows the event-by-event distribution of the 
p-r-averaged vi and V3. For X m f p /d = 2 vi has a broad 
distribution between and 0.1 with the maximum at 
0.06, which is comparable with the values at RHIC. v^s 
distribution is narrow and centered at 0.01. For smaller 
Xmfp/d ratios the viscous effect becomes smaller and the 
collective flow becomes stronger. The value of V2 and ^3 
can reach 0.2 and 0.07, respectively, for X m f p /d = 0.5, 
which corresponds to rj/s ss 0.094. Therefore, if the par- 
ton matter in high-multiplicity p+p events at LHC has 
a small rj/s ratio, both V2 and ^3 are measurable quan- 
tities. Although the Boltzmann equation is not strictly 
valid for systems with X m f p /d < 1, its solution agrees 
well with results from hydrodynamic calculations with 
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FIG. 4: (Color online) Event-by-event distribution of V2 and 
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corresponding rj/s ratio [56J. 

We are more interested in the translation from the ini- 
tial spatial asymmetry to the final collective flow. In Fig. 
[5] we show the distributions of the difference between ^ n 
and for n = 2, 3 [see the definitions in Eqs. ([3]) and 
(JBJ)]. The distributions of n = 2, 3 are almost the same. 
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FIG. 5: (Color online) Distributions (arbitrary unit) of the 
difference between the initial event-plane angle $ n and corre- 
sponding final one ty n . 

They peak at zero and have a form looking like the Dirac 
delta function. It seems that translations from £2 to vi 
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and from 63 to v 3 take place independently. 

However, the event-plane correlations present a differ- 
ent picture, which is given in Fig. [6] We have selected 
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FIG. 6: (Color online) Correlations (arbitrary unit) of initial 
£2 — £3 and final 112 — Vs event-plane angles. The dotted curve 
is obtained by independent samplings of ^2 and according 
to Fig. [5] Events with V3 > 0.005 are selected. 

events with v 3 > 0.005, because in events with v 3 < 0.005 
^3 (possibly also 4> 3 due to tiny e 3 ) is a random num- 
ber within [— 7r/3,7r/3] and has no correlation with the 
initial $3. The solid curve shows the angular correla- 
tions of the initial 62 — £3 event-planes, which is almost 
the same as the one demonstrated with a simpler initial 
condition, seen in Fig. [3] The dashed curve shows the 
corresponding correlation of the «2 — w 3 event-planes af- 
ter parton cascade simulations. Surprisingly, the V2 — w 3 
event-plane correlation function has a maximum at zero 
degree and, thus, is totally different from the 62 — e 3 
one. With the same event- by- event initial conditions but 
without parton cascade simulations we sample ^2 and 
* 3 independently according to Fig. [5] The result is 
plotted by the dotted curve in Fig. [51 which shows, as 
expected, the same trend as the initial angular correla- 
tion, although the correlation is weaker due to the width 
in the distributions seen in Fig. [5J Because the inde- 
pendent samplings of ^2 and ^3 does not reproduce the 
^2 — *3 correlation (dashed curve) in the parton cascade 
calculations, we conclude that elliptic and triangular flow 
are correlated during the dynamical expansion. This dy- 
namical correlation seems to rotate different event-planes 
(30 degree in A$) to a unified event-plane (0 degree in 
A*). 

To convince ourselves of the finding, we show the con- 
tour plot d 2 N/d(§ 3 - $ 2 )/rf(*3 - ^2) in Fig. [3 Integral 
over ^3 — ^2 gives the solid curve in Fig. |6l while in- 
tegral over $3 — $2 gives the dashed curve in Fig. [5] 
The difference between the solid and the dashed curve in 
Fig. |B] is reflected in the asymmetry in Fig. [7] along the 
plane * 3 - * 2 = $3 - $2- At a fixed $3 - <f> 2 , *3 - *2 
has a broad distribution with a center moving toward 
^3 — ^2 = 0. This is better observed in Fig. [SI where 
the final |^2 — ^3! correlations from events within sepa- 
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FIG. 7: (Color online) Contour plot d 2 N/d($> 3 - <h)M*3 
^2) (arbitrary unit). 



rate bins of the initial correlation angles, |<£>2 — $3! =0-5, 
5-10, 10-15, 15-20, 20-25, and 25-30 degree, are shown. 
We clearly see strong broadening of all the distributions 
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FIG. 8: (Color online) Correlations (arbitrary unit) of final 
V2 — V3 event-plane angles within the separate bins of the 
initial correlation angles, |$2 — $3| =0-5, 5-10, 10-15, 15-20, 
20-25, and 25-30 degree. 

toward zero degree. For instance, we choose the curve 
denoted by 10— 15 degree. When assuming independent 
translations from e n to v n , this curve is expected to peak 
at 10 — 15 degree. However, we realize its maximum at 
zero degree. 

Figure [S] shows another contour plot d 2 N/d(^>2 — 
$ 2 )/d(*3 - $3)- Integral over * 3 - $3 (or \& 2 - $2) 
gives the solid (or dashed) curve in Fig. [5] If ^2 — $2 and 
'i'a — <I>3 are independent, the contour structure should be 
symmetric along the planes ^2 — ^2 = and $3 — $3 = 0. 
It is not the case. Assuming $2 = and <f> 3 = 30 de- 
gree, (which are more favored,) and choosing a ^2 with 
^2 — $2 > 0, we see from the contour plot that events 
with ^3 — <E> 3 < are more favored. This indicates that 
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FIG. 9: (Color online) Contour plot d 2 N/d(^ 2 - $ 2 )M*3 
$3) (arbitrary unit). 



the angle between $2 and ^3 is smaller than 30 degree 
and the two event-planes rotate toward each other during 
the dynamical expansion. 

Wc have to note that our conclusion on the dynamical 
correlation between 1*2 and V3 needs further verifications, 
because the fluctuating initial configuration of hot spots 
may affect the final event-plane angular correlation. It is 
worthwhile to study this issue for a smooth initial condi- 
tion with few components e n and given initial event-plane 
angular correlation. We leave this as a task for a future 
investigation. 



IV. CONCLUSIONS 

In this work we have calculated the elliptic and tri- 
angular flow in high-multiplicity proton-proton collisions 



at the LHC energy 14 TeV. The reason for the measur- 
able flows is the assumed initial fluctuation in the hot 
spot scenario. The motivation of this work was to find 
final V2 — V3 event-plane correlation expected by the ini- 
tial 62 — £3 event-plane correlation. The latter exists for 
the assumed initial condition with three statistically dis- 
tributed hot spots originating from three independent 
fragmenting strings in p+p collisions modelled by HI- 
JING. The results obtained by using parton transport 
model BAMPS showed the largest w 2 — V3 event-plane an- 
gular correlation at zero degree, which is the opposite to 
the expectation at 30 degree. This observation indicates 
a dynamical correlation between elliptic and triangular 
flow during the expansion. Their event-planes rotate to- 
ward each other. If so, any initial correlations will be 
washed out and it is more difficult to extract initial con- 
ditions from the flow observations. On the other hand, 
because we do not expect initial event-plane angular cor- 
relations in nucleus-nucleus collisions at RHIC and LHC, 
measurements on the final harmonic flow event-plane an- 
gular correlations in these experiments would confirm our 
conclusion, if data favor zero-degree correlation of , ^3 
event-plane angles. The correlation of different flow com- 
ponents is a new finding and needs further verifications 
in the future. 
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